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Abstract. We here specialize the well known Positive Real Lemma (also 
known as the Kalman-Yakubovich-Popov Lemma) to complex matrix-valued 
rational functions, (i) generalized positive even and (ii) odd. On the way we 
characterize the (non) minimality of realization of arbitrary systems through 
(i) the corresponding state matrix and (ii) moving the poles by applying static 
output feedback. We then explore the application of static output feedback to 
both generalized positive even and to odd functions. 



1. Introduction 

For a half of a century, the Positive Real Lemma (also known as the Kalman- 
Yakubovich-Popov Lemma) has been recognized as a fundamental result in System 
Theory. We here exploit it to study two classes of rational functions, (i) generalized 
positive even and (ii) odd. 

Let C + and C_ be the open right and left halves of the complex plane respectively, 
and Pfc, (Pfc) be the set of all k x k positive definite (semidefinite) matrices. 
Recall that apx p- valued function F(s), analytic in C+ is said to be positive if 

(1.1) F(s) + F*{s) ep p seC+. 

The study of rational positive functions, denoted by V, has been motivated from 
the 1920's by (lumped) electrical networks theory, see e.g. [§],[§]■ From the 1960's 
positive functions also appeared in books on absolute stability theory, see e.g. [37] . 
[39]. 

Apx p- valued function of bounded type in C+ (i.e. a quotient of two functions 
analytic and bounded in C+) is called generalized positive QV if 

(1.2) F(iu) + F*(iuj) e P p a.e. ujeR, 
where F(ioj) denotes the non-tangential limilQ of F at the point ioj. 
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It is interesting to note that both the function set QV and its subset V, are closed 
under positive scaling, sum and inversion (when the given function has a non- 
identically vanishing determinant). We emphasize that for simplicity we adhere 
hereafter to the rational case. We shall find it convenient to use the notation 

F*{s) := F*(-s*). 

Recall that a matrix valued function F(s) is called even if 

(1.3) F(s) = F*(s). 

We shall denote by Even the set of even functions. Having F £ Even in particular 
implies that its zeroes and poles are symmetric with respect to the imaginary axis. 
Moreover, on the imaginary axis this F is Hermitian, i.e. 

F(ioj) = F*(oj) a.e. w G M. 

We shall say that apx p- valued functiorQ ^(s) is even generalized positive, denoted 
by QVE, if it satisfies both (TOJ and CCS]), i.e. 

(1.4) ePp a.e. wel. 

Namely, QVE = QVf] £ven. Note having that * £ QV is equivalent to (* + *#) G QVE . 
Scalar rational QVE functions were recently studied and then applied to Nevanlinna- 
Pick interpolation in [H Section 5]. 

In analogy to Even functions in (jl.3p . we shall say that a function F(s) is odd, 
denoted by F £ Odd, if 

(1.5) F{s) = -F*(s). 

This implies that on the imaginary axis Odd functions are skew-Hermitian, 

(1.6) F(ioj) = -F*(iu) a.e. wel. 

Recall also that a matrix valued function F(s) can always be partitioned to its 
even and odd parts, i.e. 

F(s) = F cvcn (s) + F odd {s) F even :=j(F + F#) F odd := i(F - F#) . 

Clearly, F cvcn G Even and F odd £ Odd. From (|1.2p and (| 1 .€>[) it follows that in fact 
Odd C QV. Scalar rational Odd functions were recently studied in [H Section 4]. To 
further motivate the study of Odd functions, recall that the classical Nevanlinna- 
Pick framework it was shown in |45j that if there exists an interpolating function 
within V, without loss of generality (but compromising the minimal degree) there 
exists an interpolating function within VO := Vf]Odd. For an application of this 
observation see [TU Corollary 5.2.2]. 

Specializing the above even-odd partitioning to ^ £ QV one obtains 

*( S ) = *cvo„(s) + *odd(s) *cvon G QVE * odd G Odd. 

As already mentioned, scalar, rational QVE and Odd functions were recently studied 
in [4j Section 5] and [4j Section 4], respectively. 



To ease reading, hereafter F(s) denotes an arbitrary rational function and 4 , (s) is from the 
subset of QV. 
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We now recall in the concept of state space realization. Let F(s) beapx p-valued 
rational function analytic at infinity, i.e. lim F(s) exists. Denote by q the 

s — > oo 

McMillan degree of F. Namely, F admits a state space realization 

(1.7) F(s) = C(sI-A)^B + D L:={£B) 

with A G C nxn , n>q, B,C* G C nxp and D G C pxp , namely, L G £(™+p)x(™+p) . 
If the McMillan degree of F(s) satisfies q = n, the realization is called minimal. 

The aim of this work is to characterize the realization of rational QV£ and Odd 
functions and to give an application to static output feedback. 

Related works on realization for the non-rational case exist, see e.g. [16| and more 
recently [7]. Restricting the discussion to rational functions, enables to offer simple 
explicit formulas. On realization, our main result is as follows. 

Theorem 1.1. Let *£(s) be a p x p-valued rational function so that lim M/(s) 

s — > oo 

exists and let L be the associated realization matrix, 

V(s) = C(sI-A)- l B + D L:=(g g), 

see (|1.7p . Let q denotes the McMillan degree. 

A. Let L G c(2«+p)x(2n+p) yjfth q even, 2n > q. The following are equivalent. 

(i) VP e QVE. 

(ii) There exist 2n x 2n Hermitian non-singular matrices Hi and H2 so that, 

HL + L*H g F 2n+p ML + (ML)* = 2n+p , 

where 

H := diag{£Ti, I p } M := diag{# 2 , il p }. 

(iii) There exist realization matrices L whose sub-blocks are 

(1.8) A =(i -A') B =(-d') C=(C DB*) D = DD* 

where A is n x n, B,C* are n x p and D is p x p. Furtheremore without 
loss of generality one can take spectrum A to be in C+ . 
Moreover, if the realization is minimal, i.e. 2n = q, up to similarity, L is 
of the above form. 

B. Let L G £( n +p) x ( n +p) w ith n>q and let v be so that A, the nxn part of L has 
at most v eigenvalues in C_ and at most n—v eigenvalues in C + . 

Consider the following statements. 

(i) * e Odd 

(ii) There exists a non-singular nxn Hermitian H s.t. H := diagjiJ , L p } 
satisfies 

HL + L*H = n+p . 

(iii) Up to similarity the sub-blocks of the realization matrix L are 

(1.9) A =(a4 2 ) B = (b 2 ) C={Bl-B* 2 ) D = T 3 

with T\,T 2 , T3 skew- Hermitian (Tj = —T*) of dimensions vxv , (n — v) x (n — v) 
and p x p, respectively. 
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Then (ii) and (iii) are equivalent and in turn imply (i). 

If in addition the realization is minimal, i.e. q — n, the the converse is true as well. 

The result of part A can be compared with [8] Theorem 10.2] where a state space 
realization of QVE functions is given. While the development in [8] is self contained, 
we here rely on two classical results: Lemma l2~Tl and Theorem l2 . 2l below. This allows 
us to obtain the above Theorem 11.11 A which is more general than [8j Theorem 
10.2] in the following sense, (i) minimality of the realization is not assumed and 
(ii) D — lim ^(s) can be of any rank (including zero). Moreover, the explicit 

s— >oo 

formulation in (jl.8j) and (jl.9j) enables us to directly address applications like static 
output feedback. For example, recall that a QV£ function which is analytic on iK 
admits a spectral factorization GG# with G(s) analytic in C+, see e.g. [8j Chapter 
9], [20] . [31] Section 19.3] and [3D] . We here state a result, where some of the details 
will be clarified in the sequel. 

Proposition 1.2. Let ^(s) be a p x p-valued rational QVE junction, which is not 
analytic on iM. Assume that lim ^(s) = and let the state space realization be as 

s— >oo 

in (TlU) with D = 0. 

There exists a static output feedback gain K , —K £ V p so that the closed loop 
system (I — )~ \P is analytic on iWL, if and only if for all r el the two following 
matrices 

(, .„.:„) «* ( A i In ) 

are of full rank. 

Clearly, it is suffices to check the conditions only for all ir £ spect(A). 

The outline of the paper is as follows. Section[2]is devoted to providing a perspective 
on relevant existing literature and background to be used in the sequel. Part A 
of Theorem 11.11 is proved in Section [3] Aspects of the (non) minimality of the 
realization in (jl.8[) are addressed in Section 01 On the way, we introduce a test for 
the non-minimality of a state space realization of an arbitrary system (vanishing at 
infinity) by examining common eigenvalues between the (n + p) x (n + p) system 
matrix L and its n x n submatrix A, see (jl.7p . Part B of Theorem 11.11 is proved 
in Section [5l In Section [6] we relate minimal realization of an arbitrary system 
(vanishing at infinity) with the ability of moving its poles through static output 
feedback. As a sample application of Theorem ll.il in Section [7] we study the effect 
of static output feedback on QVE and Odd systems and then prove Proposition ll.2l 

2. Background and perspective 

In this section we state known results to be used in the sequel. Generalized positive 
functions were introduced in the context of the Positive Real Lemma (PRL), see 
[5] and references thereirH. Applications of QV functions to electrical networks 
appeared in [26] , and to control in |36j , where they first casted in a Linear Matrix 
Inequality (LMI) framework, see e.g. [10] for more information on LMI. For more 
application of the generalized PRL, see [22] . 



The original formulation was real. The case we address is in fact generalized positive and 
complex, but we wish to adhere to the commonly used term: Positive Real Lemma. 
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We now recall in three characterizations of QV functions. We start with the Positive 
Real Lemma (PRL) as first presented in [T4] Theorem 1] (up to substituting the 
real setting by a complex one f| For sample of other versions the PRL see e.g. [5] 
Theorem 15.2], [22] Section 3A] and the discussions in [3] Section 2] and in [TO] p. 
34]. 

Lemma 2.1. Let L be a (n + p) x (n + p) realization matrix of a p x p-valued 
rational junction $ (s) of McMillan degree q, see (|1.7| . 

(I) IP € £/"P £/ there exists a Hermitian non-singular H 6 C™ xn so f/j a f 

(2.1) HL + L*H e¥ n+p H = Amg{H, I p }. 

If in addition the realization is minimal, i.e. q — n, the converse is true as well: 
* G QV implies (f2~T) . 

(II) If in part (I) -H G P„ then *eP. 

It is of interest to recall that in [3j Section 7] we pointed out that if L satisfies (|2.1|) , 
whenever non-singular, also satisfies (|2.1[) . with the same H . 

As a second characterizations of QV functions recall that apx p-valued function 
\l/(s) belongs to QV if and only if it can be factored as 

(2.2) *(s) = G{s)P(s)G*{s), 

where G, P are p x p-valued, G analytic in C_ and P € V . See for instance [21] for 
the corresponding result in the setting of functions meromorphic in the open unit 
disk rather than in the right open half-plane, see the discussion in [31 Section 1]. 
To present the third characterizations of QV functions we briefly mention that a 
rational function ^ is in QV if and only if the kernel 

*(s) + *(w)* 
s + w* 

has a finite number of negative squares in its domain of definition in C+. This 
equivalent characterization of (not necessarily rational) functions of the form (12.21) 
appeared for the scalar case in [13] and [15] and extended in [33], [34], [35]. The 
significance of (|2.2p to scalar rational QV functions was recently treated in [2] and 
[3] , where a more complete survey of the literature can also be found. 

The following result is crucial to our construction. 

Theorem 2.2. Let ^(s) be a matrix valued rational function. The following are 
equivalent 

(i) * € QV£ 

(ii) * G QV and in addition ^f(iu) = a.e. w G K. 

(iii) 'P(s) = G(s)G(s)# and without loss of generality G can be chosen so that 
G and G^ 1 are analytic in C_ . 

Proof The equivalence of (i) and (ii) follows from the definition of QV£ functions, 
see ([L~i]> , 

The fact that (ii) ==> (iii) is deep and was established in [44] Theorem 2]. 
The fact that (iii) => (i) is straightforward. □ 



The original formulation was real. The case we address is in fact generalized positive and 
complex, but we wish to adhere to the commonly used term: Positive Real Lemma. 



6 



D. ALPAY AND I. LEWKOWICZ 



The problem in item (iii) of finding GG# out of \P G QVE is known as spectral 
factorization whenever '5 and 'J -1 are analytic on the imaginary axis and if this 
restriction is relaxed, pseudo spectral factorization. As sample references on spectral 
factorization, see e.g. Chapter 9] [20] [31] Section 19.3] and [40]. On pseudo- 
spectral factorization, see e.g. [8] Chapter 10], [41] and the earlier work of Youla on 
which we rely, [44] . The gap between spectral and pseudo spectral factorizations 
is addressed in [8j Chapers 9, 10]. Here, Proposition 11.21 is restated and proved in 
Proposition 17.41 below. 

For completeness we mention that the spectral factorization problem has been ex- 
tended to the case where F € Even can be factored to 

F(s) = G(s) ■ diag{7, -I}-G*(s) 
see e.g. [23], [32] and [i Part VII]. 

Note that * £ QV£ if and only if one substitutes in (|2~2|) * = vj/# thus it is 
equivalent to having in (I2.2p P{s) = P for some constant positive semidefinite P. 
This conforms well with item (iii) in Theorem 12.21 

The following result from [TJ Theorem 4.1] , (see also J3H1 Proposition 2.1]) is adapted 
to our framework. 

Theorem 2.3. Let F(s) be a matrix valued rational function admitting a state 
space realization L as in (|1.7p of McMillan degree q. If 

(2.3) ML + (ML)* = n+p 
where 

M := diag{ff, il p }, 
for some Hermitian non-singular H £ C nxn , then 

(2.4) F{iu)) = F(iuj)* a.e. wel. 

Conversely, let L be a minimal realization of F(s) in (|2.4|) . i.e. q = n. Then L 
satisfies (|2.3jl . 

Note that M is not Hermitian. 

Clearly, if F satisfies l|2.4[) . so does — F and also F _1 (whenever the determinant 
is not identically zero). It is less obvious that some of this properties hold for the 
system matrix L. Namely, if L satisfies (|2.3|) with M = diag{-ff, Up}, then this 
holds for — L and whenever exist, ±L _1 satisfies (|2.3|) with M = diag{— H, il p }. 

We conclude this section with recalling in a technical result on the state space real- 
ization of a composition of a pair of rational functions (series or cascade connection 
of systems in Electrical Engineering terminology). Namely F a (z), Fp(z) are of 
compatible dimensions and F 1 (z) is obtained by 

(2.5) F 1 {z) = F a {z)F p (z). 

Assuming the state-space realization of each F a (z) and Fp(z) is known, one can 
construct a realization of the resulting Fj(z), see e.g. [27] Subsection 8.3.3], (43j 
Eq. (4.15)]. 

Observation 2.4. Given I x q and q x r valued rational functions F a (z), Fp(z), 
respectively, admitting state space realization 

F a (z) = C a (zI-A a )- l B a + D a A a £C p « xp « 
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F fs (z)^C p (zI-A fs )- 1 B p + D A^eC*^ 
B l3 eC p e xr C p e C qxp f D f3 eC qxr 

The system matrix L 1 associated with the realization of F^(z) in (|2.5p is given by, 





( A a 


B a C/3 


B a D p 








Bp 




\ c a 


D a Cf3 


D a D p 



3. Proof of Theorem 11.11 A 

(it) 

Substituting in Theorem 12.21 (it) =>■ (i) a combination of Lemma \2. II (the PRL) 
and Theorem 12.31 establishes this part. 



(i) (Hi) 

First, substitute in (fT7f|) F(s) 
tion is 

(3.1) 



G(s) so that its (not necessarily minimal) realiza- 



L„ 



A 


B 


C 


D 



Then the corresponding realization of G(s)& is given by, 



-A* 


-C* " 


B* 


D* 



Now, substituting in Observation ^. 41 F n = G and Fp — G# along with Proposition 
2J2I establishes the structure in (11.81). i.e. 







( A 


BB* 


BD* 


(3.2) 


-1 





-A* 


-C* 






V c 


DB* 


DD* 



Finally, we have shown above how to construct a realization of QVE function of 
dcgrcc= 2n where n is arbitrary, see (|3.2j) . Note now that if "J € QV£ is realized 

by L = 





A 


B N 


■(■ 


C 


D , 


A 


* 


B 





* 





C 




D t 



of the form in (|3.2p . then the same ^ can also realized by 
, where * means an arbitrary block of suitable dimensions. It 



may be the case that one can not transform L by change of coordinates to the form 
in (|3.2[) (in particular the dimension of the A part in L may be odd). However, as 
all minimal realizations are similar, see e.g. [HI 8. 29], [37J Theorem 2.4-7], one of 
them is of the form (13.21). 



(Hi) ==> (ii) 

We find it convenient to introduce the following four intermediate conditions: 



(3-3) _ 

Condition (ii) in Theorem ll.il A is satisfied where Hi and H2 share the same inertia. 
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(3.4) Condition (ii) in Theorem ll.il A is satisfied with H 2 unitarily similar to H\. 



(3-5) _ 

Condition (ii) in Theorem ll.il A is satisfied with H2 unitarily similar to H\, both involutions. 
(3.6) 

Condition (ii) in Theorem ll.il A is msatisfied with Hi = ( j 7 q ) H2 = i ( j n ) . 

Trivially ((33]) ([33]) (23) ([33J (ii), so all we need to show 

is that (iii) ([376]) . 

Indeed, a straightforward calculation reveals that the matrix L in (|3-2[) satisfies 

the conditions in (l3~6l) where iJL + L* H = ( Vzb ) ( V2B ) . Verifying (|2~3l) is 

_ \V2dJ \V2dJ 

immediate so the proof is complete. □ 



Note that the matrix L in 
(3.7) 



has a special symmetry, i.e. 





--fn 









!)( 










A BB* BD* 
-A* -C* 
C DB* DD* 



A* C* 
A BB" BD* 
C DB* DD* 



In particular, the A matrix in (II. 8|) (the upper left bock in L) has a Hamiltonian 
structural, see also (|7.2p below. This will be exploited in the proof Proposition 17.41 
below. 

We now examine an aspect of Theorem 11.11 A: Above we have indirectly proved 
that A (ii) =*> (3~3]) =>■ (E3) =>■ IJ33JI (EU) => A (iii). 

We now illustrate the fact that directly showing these implications is not obvious. 



where Pi , P 2 G P„ and P 3 G P; are arbitrary. Then A; 



(1331) =^> 


(iii). 










P2 





Let A, = 


(5 















P3 



satisfies the conditions of (|3.6[) with the same -ff and M as in (|3.6p . This in turn 
implies that also L + A/ satisfies the same conditions. However, L + A; is no longer 
of the form of (|1.8[) . This is illustrated by the following example. 

1 

1 ] of the of the form of (|1.8p be a realization the rational 
1. Now in A/ substitute Pi = 1, P2 = 2 and 
. Indeed, A; satisfies condition Q3.6P and so does 



Let L x = 1 

V 1 1 

QV£ function ipi(s) 



1 







2 





P 3 = 3, i.e. A, = 






















( -\ 


3 




L 2 := L x + Ai = 




1 








1 


4 



-1 . Although L2 is not of the form of ()1 .8[) . it is 



a state space realization of the rational QVS function i\)%(s) = j^i + 4. Indeed, 



^sometimes called H-skew-Hermitian, see e.g. 1311 sub-section 3.2]. 



REALIZATION OF EVEN GENERALIZED POSITIVE AND ODD RATIONAL FUNCTIONS 9 



ip2 (s) can also be realized by L' 2 

of CEU). 



-2 4 


4 v 


2 


2 - y/5 


VE- 2 4 


4 , 



which is of the form 



Consider the condition for the submatrix that HiA+A* Hi £ P„ and -f^A + A*H2 = r , 
Take A = ( J g). i?i = diag{l, -5} and iJ 2 = ( ^ ) • Then the conditions are 
satisfied with Hi and H2 unitarily similar, but neither is an involution. 



Consider the condition for the submatrix that Hi A+A*Hi £ P„ and H2A + A* H2 = 0„. 
Take A = (? jj). #1 = diag{l, -5} and iJ 2 = (1-3). Then the conditions are 
satisfied with iii and H2 sharing the same inertia, but they are not unitarily similar. 



4. MINIMALITY OF THE REALIZATION 

We now address the question of minimality of the obtained realization. We first 
resort to the following observation which goes beyond the scope of this work. 

Lemma 4.1. Let L £ (£(«+p) x (™+p) ^ e realization of apxp-valued rational function 
F(s) of McMillan degree n. Namely, 

F(s)=C(sI-A)- 1 B + D 

see (|1.7I) . If the realization is not minimal, there exists A £ spect(j4) so that 
A £ spect(L) for all D. 

If rank(_B) = rank(C) or the number of Jordan blocks in A is less or equal to 
min(rank(i?), rank(C)) then the converse is true as well. 

Proof We first show that if the realization is not minimal there exists A £ spect(A) 
so that A £ spect(L) for all D. 

If the realization L is not controllable, there exists 7^ v £ C" so that v* (A — XI n ) = 0, 
for some A £ C and in addition v* B = 0, see e.g. (BJ Theorem 3.3.1], [27l Theorem 2.4-8 
[3T1 Theorem 4.3.3]. This implies that (n + p)-dimensional vector v : =(„) satisfies, 
v*(L — XI n+p ) = with the same A. Namely, A £ spect(A) f] spect(i). 

Similarly, if the realization L is not observable, there exists 7^ v £ C" so that 
(A — \I„)v = 0, for some A £ C and in addition Ci — 0, see e.g. [U Theorem 3.3.6], 
[27l Theorem 2.4-8]. This implies that (n + p)-dimensional vector v :=uj satisfies, 
[L — XI n+p )v = with the same A. 

For the converse direction we now show that if the realization L is minimal one can 
construct D so that the matrix 

T \T _ ( A-\I n B \ 
Li Al n+p — I q D-XIp J 

is nonsingular for each A £ spect(A), i.e spect(A) f] spect(L) = 0. To this end, take 
D so that spect(A) |~| spect(Z?) = 0. Hence, one can write 

r xr —fin B(D-\I p )- 1 \ ( A-\I n +B{XI p -D)- 1 C \ 
L M n+p ~ [ o Ip ) { c D-XI P j ■ 
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In this case A £ spect(L) if and only if the matrix 

A - XI n + BKC with K := (D - XIp)- 1 

is singular. The construction for a minimal realization of such K £ <C pxp so that 
A — XI n + BKC is nonsingular for all A € spect(yl) is established in Proposition [(0] 
below, so the proof is complete. □ 

Admittedly, we do not know whether or not the conditions on rank(f?) and rank(C) 
in the above result are inherent to the problem or just a by-product of the technique 
we employed. 

We now illustrate the fact that it may be the case that spect(A) f] spect(L) 7^ 
for some D, although a realization is minimal. Indeed, take A — (scalar) 
B* , C £ C 2 , both non-zero, and D £ C 2x2 . Namely the following 3x3 system 
matrix L — is a minimal (degree=l) realization of a 2 x 2-valued rational 

function F(s) = -CB + D. Obviously, minimality is independent of D. However, if 
for example F(s) = (- + y)CB, namely D = 7CS, where 7 £ C is arbitrary, both 
A and L are singular, i.e. have a common eigenvalue (=zero). 

The special structure of D in the above example suggests that if a realization is 
minimal, the matrices L and A typically do not share common eigenvalues. 

One can now exploit the special structure of the realization QV£ functions in (|3.2[) . 
(|3.7j) . in applying Lemma |4~T1 to obtain the following. 

Observation 4.2. Let L g (s) in (|3.1j) be a realization of p x p-valued rational func- 
tion G(s) and without loss of generality assume that spect(^4) C C+. 
Let G(s) — G(s) — D, i.e. G(s) = C(sL n — A)~ 1 B. The following is true. 

A. Minimality of the realizations of G(s) and of G{s) is equivalent. 

B. The following are equivalent. 

(i) The realization of G{s) is not minimal. 

(ii) There exists X € spect(A) so that for all D, X £ spect(L g ) in (|3.1[) . 

(iii) For all D the realization L in (|3.2[) . of G(s)G(s)# is not minimal. 

(iv) There exists X £ spect(A) 1J spect(— A*) so that for all D, X £ spect(i) in 

(E2J. 

C. Assume that the realization of G(s) is minimal. If spect(A) (£ iM., one can 
always find D so that the realization of G(s)G(s)^ in (|3.2|) . is not minimal. 

Proof A. As minimality of realization is equivalent to controllability and observ- 
ability, the claim is obvious. 

B. First note that if C(sl — A)~ 1 B + D is a realization of a QV£ function of the 
form ([321), then from §ZJ\, it in particular follows that C* =(£ ~^)B, which 
implies rank(B) = rank(C). Thus, one can apply Lemma 14.11 which now implies 
(i) <^=> (ii) and (iii) <^=> (iv). 

(i) => (iii) 

If the realization of G(s) is not minimal, assume it is not observable. Namely, 
there exists O^vG C" so that ( * ) is an (n + p)-dimensional vector of Lq in (13. ip 
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(corresponding to an eigenvalue A). Then, independent of D, is an (2n + p)- 

dimensional vector of L in (|3.2|) (corresponding to an eigenvalue A). Thus, this 
realization is not observable, (controllability can be similarly treated). 

(iii) => (i) 

Assume that the realization L in (13.21) is not minimal for all D. If it is not observable 



(controllability can be similarly treated), there exist v±,V2 € C" so that \ vlj is a 

non-zero (2n + p)-dimensional eigenvector of L, corresponding to an eigenvalues A. 
As this holds for all D, it implies that B*V2 — 0. Now if V2 7^ then it is also 
a left eigenvector of A (corresponding to an eigenvalues —A*) so the pair A,B is 
not controllable. If = then 7^ v\ is an eigenvector of A (corresponding to an 
eigenvalues A) and Cv\ — 0, thus the pair A, C is not observable. 

C. Let A 6 C+ be an eigenvalue of A and let 0/»€ C™ be a corresponding left 
eigenvector. Namely, v*(A - XI) = 0, i.e. (A* — X*I)v — 0. Minimality of the 
realization of G(s) implies that v*B 7^ and Cv 7^ 0. 

By construction the matrix (A + X*I) is nonsingular, so the non-zero (2n + p)- 
dimcnsional vector ( ( " + J j is well defined. Now, for D — C(A+ X* I)^ 1 B 

this is an eigenvector of L in (|3.2p (corresponding to an eigenvalue —A*) and thus 
this realization is not observable (controllability can be similarly treated). □ 

We conclude this section by illustrating the fact that the condition spect(^4) <f_ iM. 
in item C is essential. 

Example 4.3. If spect(A) C iR, it may be that the realization L in (|3.2ll is minimal 
for all D. Take for example the scalar function 'I'i(s) = - which is positive and 
odd. Its minimal state space realization is ^i(s) = C(sl — A)^ 1 !} with A = 0, 
5 = 1 and £7 = 1. 

Let * 2 € QVE be of the form * 2 (s) = *i(s)*f (s) = —js- Its realization is of the 
form of (|1.8[) . given by 

(4.1) A=(gS) S = (_1) £7 =(1,0) D = 0. 

In particular, this realization is minimal. 

Consider now the QVE function ^(s) + D, where D > is a parameter. For all 

D, the realization is minimal, i.e. of degree 2. □ 



5. Odd systems 

Proof of Theorem [□] B 

(ii) => (iii) 

Let has ^ eigenvalues in C_ and n — v eigenvalues in C+ . Then up to similarity 
on L (and A) and congruence on H (and H) one can take H = diag{— I v , I n - V } 
(and thus H = diag{— I v , I n -v, Ip})- Substituting in HL + L* H = 0„ +p yields 
CU). 

(iii) => (i) 
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A straightforward calculation yields ^(s) = C(sl n — A) 1 B + D. We shall concen- 
trate on ^(s)| s=iiij for We first conformally partition 



(it si A)- 1 - ( ^^- T i -A \ 1 ._ fW(iu) X(iu)\ 
(lUI n -A) _ A , iulIn u _ T2 ) — [ y(iu) z(iu ) ) ■ 

We thus formally have, 

(5.1) = B*W{iu))B x + B\X(iu)B 2 - B'*Y(iuj)B 1 - B* 2 Z{iu)B 2 + T 3 . 

As for all us E M the matrix i ■ diag{/„ , — I n _ v }(iu;I n — A) is Hermitian, so is its 
inverse 



(iwl n - A) 1 • i ■ diag{-/„ , I n - V ] 



-iW(iui) iX(iui) 
— iY(iuj) %Z{iuj) 



Namely, W(iu) = -W(iu))* := T 4 (iuj), Z{iu) = -Z(iu)* := T 5 (ioj) and X(iuj) = Y(iu)* 
Substituting now in (|5.ip yields 

= S;r 4 (tw)Bi + B*X(iu)B 2 - (BlX(iu))B 2 )* - B* 2 T 5 (iuj)B 2 + T 3 . 

Namely (|1.6j) is obtained and this part of the claim is established. 



(i) => (ii) 

First consider the condition in (ii) and note that if one denotes N := diag{/„, il p } 
then 

(5.2) = HL + L*H = HNN*L + L*NN*H = ML + L*M* = ML + (ML)*, 

where M := diag{7J, il p } and L = N*L =(f c _f D j i.e. if L realizes 

\fr(s) = C(sl n - Ay 1 B + D, L realizes F(s) := -i$(s). It now follows that on 
the imaginary axis F(s) is Hermitian, see (|2.4j) and ^(s) is skew-Hcrmitian, i.e. 
W e Odd. 



Assuming the realization is minimal, i.e. q = n, both directions of Theorem 12.31 
holds, so in the proof is complete. □ 

Note that as Odd C QV, Lemma [2~Tl (the PRL) is satisfied: Indeed f|5.2[) is a special 
case of (HID, 



Of a special interest is the set VO of positive-odd functions (a.k.a Lossless or Fos- 
ter), i.e. Oddf]P. In electrical networks theory they are associated with L-C cir- 
cuits. For more details see e.g. [9j 2.20, 2.36, 2.39, 7.33, 8.35, 8.36, 8.37], pp. 
12, 221, Theorem 2.7.4], Q21 subsections 4.2, 5.1]. Recall also from ([221) ^ follows 
that * e Odd can always be factored as * = GPG* with P G PC Thus in a 
sense, PC functions generate all Odd functions. The following is known, see e.g. 
[51 Eq. (5.2.6)]. It is immediate from part B of Theorem ll.il upon substituting in 
PI) v = n. 



Corollary 5.1. Let^(s) be a pxp-valued rational function so that lim \P(s) exists. 

s— *-oo 

Let L £ (£(™+p)x(™+p) 5 e £/j e associated system matrix. If 
L = 







T„ec" x " t„=-t* 


B* 


T P ) 


T p ec pX! > t p =-t; 



T/ierc, * € PC 

J/ i/ie realization is minimal, up to similarity, the converse is true as well. 
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We conclude this section with the following observation. Clearly if vE'i(s), ^(s) 
are VO functions of the same dimensions, then ("Ji — ^2) € Odd. Comparison 
between the above Corollary and part B of Theorem 11.11 reveals that not every 
Odd function can be written as a difference of VO of functions. 



6. Static output feedback - arbitrary functions 

Recall that applying a static output feedback to an input-output system y{s) — F(s)u(s) 
(it is m-dimensional input and y is p-dimensional output) means taking u = Ky + u' 
with it' an auxiliary input and K a m x p constant matrix. The resulting closed 
loop system is y(s) = F c .i.(s)u'(s) with F c ,\. = (I p — FK)~ 1 F. 

For simplicity, we adopt the common assumption that lim F(s) = ("strictly 

s— foo 

proper" in engineering jargon). Thus, in (|1.7j) the (n + p) x (n + m) realization 
matrix L is of the form L = ( ^ | ^ ) • After applying a static output feedback, the 
closed loop realization matrix L c \ is 

(6.1) Lcl. = ( A g B Q ] U := A + BKC. 

The simplicity of the static output feedback has made it very attractive. However, 
exploring its properties turned out to be challenging. The most common associated 
problem has been stabilization, namely guaranteeing (in the continuous time case) 
that spect(A c i) C C_. This is illustrated in a basic way in [27 , Section 3.1] and for 
a sample of more recent references see e.g. [TH], [53], [3D], [35], [32]. One can go 
beyond stability, and in [3] Proposition 8.1] we characterized systems which may 
turned to be QV, and in particular V, through static output feedback. 

Our first result here goes beyond the scope of this work. To this end, we need the 
following notation. Let L be realization of a p x m-valued rational function F(s) 
of McMillan degree n with lim F(s) — 0. Namely, 

s— ► 00 

(6.2) F(s) = C{sl - A)^ 1 B L: 

see (|1.7I) . I shall find it convenient to denote j3 := rank(_B), 7 :— rank(C). 
Namely, there exist nonsingular matrices Rt € i[2 mxm and R c E C pxp so that 
R c F(s)R b =(^i s) n ) where F(s) is 7 x 8- valued. 
We shall also denote by r the number of Jordan blocks in A. 

Proposition 6.1. Let L be realization of F(s) as in (|6.2p . If this realization is not 
minimal, there exists A G spect(A) so that for all K G C mXp , X G spect(^4 c i), see 

(EH). 

If min(/3, 7) > r or 8 = 7, then the converse is true as well. 

Proof We first show that if the realization is not minimal then there exists 
A G spect(^4) so that for all K, A G spect(A c i). If the realization is not control- 
lable, there exists a left eigenvector 7^ vl G C™ (the subscript stands for "left") 
so that v* L {A — XI n ) = 0, for some A G C and in addition v* L B = 0, see e.g. 
[27] Theorem 2.4-8], [3TJ Theorem 4.3.3]. This implies that 

v* L (A cl - A7„) = v* L (A + BKC - A/„) = v* L (A - \I n ) = 0. 



A 


B ' 


C 
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Similarly, if the realization is not observable, there exists 0/i)ji£ C™ (the sub- 
script stands for "right") so that (A — XI n )vR = 0, for some A e C and in addition 
Cv R = 0, see e.g. [UJ Theorem 2.4-8]. This implies that 

{A cl - XI n )v R = {A + BKC - XI n )v R = (A - XI n )v R = 0. 

For the converse direction, assuming the realization is minimal, we now construct 
K e C mxp so that spect(A) f| spect(A d ) = 0. To this end, K should be so that 

0^»t6C" v* L (A-\I n )=0 =^ v" L (A cl -XI n )=v* L BKC^O 
0^v R eC n (A-M rl )v R =Q =^ (A ci -\I n )v R =BKCv R ^0 

where the subscripts L and R stand for "left" and "right". We shall find it con- 
venient to denote K := 5K where ^ 5 e C will be later determined. Thus, we 
actually look for K (which will turn to be an isometry if m > p, else coisometry) 
so that 

(6-3) 

O^wkSC" (A-\I n )v R =Q BKCv R ^0. 

Consider now the singular values decomposition, see e.g. [Ml Theorem 7.3.5], of 
the matrices B and C 

B = U b ^ b W b C = U C ^ C W C 

with U h , W c G C" xn , W b G C mxm and U c G C pxp all unitary and 

with E(, and E c positive diagonal of dimensions /3 x ft and 7x7 respectively. 
If /3 = 7 then it is sufficient to take 

* = w ?(oo (m _„V»)^ 

so taht 

BKC = 5U b ( tb ^ n W 6 . 

Minimality guarantees that (|6.3[) is satished. 

Next assuming that min(/3, 7) > r let 

Vl:={S*«z, : v* L (A-XI n )=0} V R := {Cv R : (A - XI n )v R = 0}. 

Namely, Vl C C m and V R C C p can be spanned by r vectors and let Jfi, G C mxr 
and JTr G C pxr be isometries whose range is equal to the span of Vl and V R , 
respectively. Thus, if one takes 

K = K L K* R , 

minimality guarantees that (|6.3j) is satisfied. Thus, there is only the value of the 
scalar S yet to consider. 

If A = XI n for some A G C, minimality of the realization implies that m > (3 = n 
and p > 7 = n. Thus, K — W b *U* so any non-zero S will do. Hence, assume 
hereafter A 7^ XI n . 
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In fact, we claim that for any \S\ sufficiently small, K satisfies the requirements. 
To this end recall, see e.g. [211 Corollary 7.3.8], that if M and A are two matrices 
of the same dimensions then for all k, 

||A|| 2 > \a k (M + A) - a k (M)\ 

where o~k{-) are the respective singular values (<7i > a 2 > " " )• For j = 1,2, ... let 
us denote by rjj the smallest positive singular value of A — Xjl n where Xj € spect(A). 
Next let rj :— min(r/i , 772 , . . . ). Taking, 

^=(2r / ||S|| 2 ||C|| 2 )- 1 

guarantees that none of the positive singular values of A — Xjl n with j = 1, 2, 
was moved to zero. Namely, A c \ — XI n is nonsingular for all A £ spect(A), so the 
construction is complete. □ 

As already mentioned, we do not know whether or not the above conditions on (3, 7 
are inherent to the problem or just a by-product of the technique we employed. 



7. Static output feedback - QVE functions 

We start with the following question: Under what conditions the QV class, and its 
subsets of QVE and Odd, are invariant under static output feedback. 

Proposition 7.1. Consider the system y(s) = ^(s)u(s) where ^ G QV is p x p- 
valued rational function s.t. lim ^(s) = 0. Consider the static output feedback 

s— > 00 

u = Ky + u' where K is a constant p x p matrix with u' auxiliary input. 

A. for all K so that —(K + K*) £ P p also the resulting closed loop system is in 
QV. 

B. If in addition ^> is even, i.e. %j) 6 QVE then, for all K so that —K £ P p also the 
resulting closed loop system is in QVE. 

C. If ^ £ Odd, then, for all K so that K + K* = also the resulting closed loop 
system is in Odd. 

Proof Item A is part of [31 Proposition 8.1 (hi)]. 

B. Substituting (|1.8[) in (|6.1I) yields, for the open loop 

(7.1) A=(^|!) S=(J.) C=(C 0) D = 0. 
Thus, for the closed loop 

Hence, for all —K £ P p , this L c .i. satisfies the conditions in (|3.6j) . so this part of 
the claim is established. 

C. Substituting {L9} in dHJ) yields, 

Thus, for all K skew-Hermitian (K = —K*) this L c j, satisfies condition (ii) in 
Theorem II. II B with H = diag{— I v , I n - V , I p }, so the claim is established □ 
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Example 7.2. We here illustrate item B in Proposition 17.11 Consider again the 
QVS function ^(s) = from Example 14.31 Applying to ^(s) a static out- 
put feedback yields 1 &2c.i.( s ) = s *~+k so ^c.i. € QVS for all > fc. Indeed, 
*2c.l.(w) = ^p—j: for all wel. □ 

We next use Proposition 17.11 B to specialize the result of Proposition 16.11 to the 
class of QV systems. 

Corollary 7.3. Let \& G QV be apxp-valued rational function s.t. lim ^(s) = 0. 
Its realization matrix L = satisfies (|2.1[) and £/ie closed loop is as in (|6.1[) . 

T/iis realization is not minimal, if and only if there exists A G spect(A) so t/iat 
A G spect(A cl ), see (EI]), /or a/Z if - (if + K* ) G P p . 

Indeed, as already pointed out in the proof of part B of Observation ^. 21 from fl3.7[) 
it follows that j3 = rank(B) is equal to 7 = rank(C). Thus, both directions of 
Proposition 16.11 hold. 

Recall that in Theorem 12.21 and the proceeding discussion we pointed out that 
v & G QVS always admits pseudo spectral factorization, i.e. ^ = GG# with G(s) 
analytic in C_. It is desired for "J to admit spectral factorization, i.e. where the 
factor G(s) is analytic in C-. Recall also that spectral factorization has various 
applications in control and filtering theories, see e.g. [6], [8j Section 16.3], [El Sec- 
tion 6] , [20] , [22] , [28] , [31] , [40] . We now address the following question: Assuming 
W G QVS, does not admit spectral factorization, under what conditions can one 
apply a static output feedback so that — (I p — ^K^ 1 ^ does admit spectral 
factorization. 

Specifically, in Proposition 17.11 and Corollary 17.31 we examined the use of static 
output feedback for moving the poles of \& G QV while retaining in the QV 
class (and QVS in particular). We know that a QVS function admits spectral 
factorization if it is analytic on the imaginary axis. Thus, the problem at hand 
is actually about moving by static output feedback poles of a QVS function away 
from the imaginary axis. 

To gain intuition we first look at the scalar case. Let ^ G QVS be written as 
"J" = j) with N(s), D{s) polynomials so that, > (including +00) for all 

Assuming that N(s) and D{s) have no common roots on iR, this in fact is 
equivalent to having for all lj G R: N(iui) > 0, D(iu) > and D(ioj) + N(iuj) > 0. 
This in turn is equivalent to having, for all cj G R: N(iu>) > 0, D(iu>) > and 
D(iu) + aN(iuj) > for all a > 0. Next recall that I'd = (I p - ^k)- 1 ^ = D N kN - 
In conform with part B of Proposition [7T] for all — k > 0: For all w G R the closed 
loop denominator satisfies (D(iuj) — kN(i<jj)) > 0. Thus in fact a — —k and the 
closed loop denominator is positive on iR. To conclude, if N(s) and D(s) have no 
common imaginary roots0, for all —k > 0, is analytic on ?R. 

Before stating the result we need some preliminaries. For a pair A G C" xn , B G C" XI 
and a region fi C C, one can define the condition, 

rank(,4 — Ai„ : B) = n V A G O. 



^In the spirit of Proposition 16.11 when restricted to the imaginary axis, the realization is 
minimal. 
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If Q = C the pair A, B is said to be controllable. If ft is a subset of C, this is 
equivalent to the existence of a matrix R £ C mx " so that spect(^4 + BR) f]Cl = 
(in particular, for ft = C+ the pair is said to be stabilizable, see e.g. [27l p. 205] 
pTTl sub-section 4.4]). 

Similarly, for a pair A £ C nx ", C £ C?> x ™ ne can define the condition, 
rank( A ^ / ")=n V A e ft. 

For ft = C the pair A, C is said to be observable, for ft = C+ the pair is detectable, 
see e.g. [311 sub-section 4.4]). 

It is of interest to point here out that in Proposition 16.11 we required minimality 
of realization, i.e. ft = C. From [30l Theorem] it follows that having the pair 
A, B stabilizable and the pair A, C detectable (i.e. ft = C+) is necessarjlj for 
stabilizability of a system by static output feedback. Below, we can be "modest" 
by resorting to ft = £R and reformulate Proposition II. 21 

Proposition 7.4. Let *5f(s) be a p x p-valued rational QVE function, which does 
not admit spectral factorization. Assume that lim ^(s) = and let the state space 

s— >oo 

realization be as in (|7.ip . 

There exists a static output feedback gain K , ~K £ ¥ p so that 'J'ci(s) admits spectral 
factorization, if and only if for r £ R the two following matrices 

are of full rank. 

Proof First, recall that ^c^s) admits spectral factorization if and only if in the 
corresponding state space realization, the spectrum of A c x avoids the imaginary 
axis. Following Proposition 17.11 this in turn is equivalent to finding conditions 
on A, B and C so that there exists —K £ F p so that Ad, in (|7.2|) . will have no 
eigenvalues on the imaginary axis, i.e. the matrix 

( 7 - 3 ) Au. ~ irhn - ( A ^ _ ( /t;„ r ) Vr £ K 

is nonsingular. Note now that the nonsingularity of A c .\. —irl2n in (|7.3|) is equivalent 
to that of 

Namely, we search for —K £ P p so that for all r £ R the Hermitian matrix M is 
nonsingular. 

Next consider the Lyapunov equation Q = HM + M*H with H — diag{— /„ , /„} 
and 

Q = 2di&g{-C*KC, BB*}. 
Now for all —K £ F p , indeed Q £ P2n- This implies that the matrix M has at most 
n eigenvalues in each open half plane, see e.g. [25l Lemma 2.4.5]. Furthermore, from 
the Generalized Inertia Theorem for the Lyapunov equation, see e.g. 25, Theorems 
2.4.7, 2.4.10], it follows that there are exactly n eigenvalues in each open half plane, 
if and only if there exists K, —K £ P p so that the pair M, Q is observable. 



7 An additional Riccati type condition makes it also sufficient. 
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Recall now that the pair M, Q is not observable, if and only if, there exists a (right) 
eigenvector of M which lies in the null-space of Q. Namely, there exists 0/j)€ C 2n 
so that Mv = Xv for some A 6 R and Qv = 0. This in turn is equivalent to having, 
(M + jQ)v = Xv and Qv = 0. Note that 

2 "* \A — irl„ ) 

To summarize: The pair M, Q is observable if and only if the pair M, Q is observable 
and this is equivalent to the nonsingularity of M (which in turn is equivalent to the 
nonsingularity of A c .i. — irlzn)- 

Now, the pair M, Q is observable if and only if there exists K, —K G P p so that 

for all r G E, both matrices: ( A ~t rIn ~ ) and ( ^ A ~} r J n) * ) are of full rank. 

\-c*kcJ \ bb* J 

Next, the matrix ^ ^ is of full rank, if and only if the matrix — irl n : B 

is of full rank, so the first part of the condition in the claim is established. 

Note now that if the matrix ^ ^3*^3 ) * s °^ ran ^ f° r some ~K G Pp it is 
of full rank for K = i.e. the matrix ( A Z tr I n ) is of full rank. This in turn is 

equivalent to having the matrix ^ ^-^ r/ » J of full rank, so the second part of the 
condition in the claim is established and the proof is complete. □ 

We now illustrate the result of Proposition 17.41 

Example 7.5. Consider the QVE function ^(s) = ^ from Examples 14.31 and !7. 21 
As shown in Example 14.31 ^2 admits only pseudo-spectral factorization. However, 
^20.1.(5) = ip-fk w ^h > ^7 see Example 17.21 can be factored to ^% c .\. = G2GJ 
with C?2(s) = -j=frr s ■ From Example 14.31 it follows that 

A «=Uic *S) = (-l 1 o) 0>k. 



Indeed, the realization A, B, C given in Example 14.31 is minimal and hence in 
particular the conditions in Proposition 17.41 are satisfied. Thus, indeed the matrix 
A c \ has no imaginary eigenvalues. □ 



From the proof of Proposition EH it follows that if <i> G QVE and K, -K G P p are 
so that = (I p — fyK)~ lx & admits spectral factorization, the same is true for aK 
where a > may be arbitrarily small. Namely, ^c^s) may be a small perturbation 
of*(s). 

We conclude this section by noting that upon comparing Propositions 16.11 and 17.41 
one can make the following statement. 

Corollary 7.6. Let G(s) = C(sl n — A)^ 1 )B be a realization of a rational function 
so that both B,C G <C pxn are of the same rank. The following are equivalent. 

• There exists a static output feedback gain K G C pxp so that the closed loop 
system (I p — GK)~ 1 G is analytic on iE. 

• There exists a static output feedback with —K G P p so that the closed loop 
system (I p — GG&K)~ 1 GG& is analytic on iE. 
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